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Abstract. Virtual corrections for electron-positron annihilation into one real and one off-shell photon of 
invariant mass are evaluated. Special attention is paid to those configurations where the real photon is 
collinear with the beam direction. This calculation is an important ingredient of a Monte Carlo program, 
which simulates events with tagged photons from initial-state radiation, including NLO corrections. 



1 introduction 

The total cross section for electron-positron annihilation 
into hadrons is one of the fundamental observables in par- 
ticle physics. Its high energy behaviour provides one of the 
first and still most convincing arguments for the point- 
like nature of quarks. Its normalization was evidence for 
the existence of quarks of three different colours, and the 
recent, precise measurements even allow for an excellent 
determination of the strong coupling at very high and 
intermediate energies (e.g. and refs. therein) through 
the influence of QCD corrections. 

Weighted integrals over the cross section with prop- 
erly chosen kernels are, furthermore, a decisive input for 
electroweak precision tests. This applies, for example, to 
the electromagnetic coupling at higher energies or to the 
anomalous magnetic moment of the muon. 

Of particular importance for these two applications is 
the low energy region, say from threshold up to centre-of- 
mass (cms) energies of approximately 3 GeV and 10 GeV, 
respectively. Recent measurements based on energy scans 
between 2 and 5 GeV have improved the accuracy in part 
of this range. However, similar, or even further improve- 
ments below 2 GeV would be highly welcome. The region 
between 1.4 GeV and 2 GeV, in particular, is poorly stud- 
ied and no collider will cover this region in the near future. 
Improvements or even an independent cross-check of the 
precise measurements of the pion form factor in the low en- 
ergy region by the CMD2 and DM2 collaborations would 
be extremely useful, since this dominates in the analysis 
of the muon anomalous magnetic moment. 



Experiments at present electron-positron colliders op- 
erate mostly at fixed energies, albeit with enormous lumi- 
nosity, with BaBar and BELLE at 10.6 GeV and KLOE 
at 1.02 GeV as most prominent examples. 

This peculiar feature allows the use of the radiative 
return, i.e. the reaction 



"'^(Pi) + e (P2) — > 7(fc) + 7*(Q)(^ hadrons) 



(1) 
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to explore a wide range of in a single experiment p,H, 

Nominally masses of the hadronic system between 2m,r 
and the cms energy of the experiment are accessible. In 
practice it is useful to consider only events with a hard 
photon — tagged or untagged — to clearly identify the 
reaction, which lowers the energy significantly. 

The study of events with photons emitted under both 
large and small angles, and thus at a significantly en- 
hanced rate, is particularly attractive for the tt+tt^ final 
state with its clear signature, an investigation performed 
at present at DAf^NE [|o[|y,|l|l. Events with a tagged 
photon, emitted under a large angle with respect to the 
beam, have a clear signature and are thus particularly 
suited to the analysis of hadronic final states of higher 
multiplicity |p^ . 

The inclusion of radiative corrections is essential for 
the precise extraction of the cross section, which is neces- 
sarily based on a Monte Carlo simulation. A first program, 
called EVA, was constructed some time ago It includes 
initial-state radiation (ISR), final-state radiation (FSR), 
their interference, and the dominant radiative corrections 
from additional collinear radiation through structure func- 
tion techniques |Q . 

The complete NLO corrections have recently been im- 
plemented in a program called PHOKHARA ||l5|. Both 
programs, however, were designed to simulate reactions 
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with tagged photons, i.e. at least one photon was required 
to be emitted under large angles. 

An important ingredient in the extension of the NLO 
Monte Carlo program PHOKHARA to small photon an- 
gles is the evaluation of the virtual corrections to the re- 
action (|l|) in the limit m\/ s <^ 1, which are equally valid 
for large and small angles. Compact results for the one- 
loop two-, three- and four-point functions that enter this 
calculation can be found in the literature for arbi- 

trary values of m^/s. However, the combination of these 
analytical expressions with the relevant coefficients is nu- 
merically unstable in the limit of small mass and angles. A 
compact, numerically stable result, valid for an arbitrar- 
ily small photon angle, is therefore required. As a conse- 
quence of the highly singular kinematic coefficients, terms 
proportional to and even m"^ must be kept in the ex- 
pansion, which after angular integration will contribute to 
the total cross section even in the limit ml/s 0. 

The present paper extends the analysis of Ref. |Q 
where the corrections from virtual and soft photons were 
presented for the case of large angles. In section 2 we re- 
call the basic definitions and describe the systematic pro- 
cedure used in the expansion of the results for small ml / s 
and small angles simultaneously. The analytic results for 
real and imaginary parts of the leptonic tensor, expressed 
in an angular momentum basis, are presented in section 3 
and compared with results for related quantities that can 
be found in the literature. After summation over the po- 
larizations of the virtual photons, our result agrees with 
the one of Berends, Burgers and van Neerven j2^,^. The 
result of Kuraev, Merenkov and Fadin |^ for the real 
part of the tensor, which was obtained for virtual Comp- 
ton scattering 7* + ^ 7 + is related to our case 
by proper analytic continuation; and indeed after analytic 
continuation we find agreement for the real partQ. Section 
4 contains our summary and the conclusions. The mass- 
dependent terms proportional to and ml, expressed in 
the Cartesian basis, are given in Appendix^ The scalar 
loop integrals needed for the calculation are listed in Ap- 
pendix H. 



2 The leptonic tensor for the radiative return 

Consider the e+e^ annihilation process 

e+bi) + e-b2)^7*(Q)+7(fci) , (2) 

where the virtual photon decays into a hadronic final state, 
7*(Q) —* hadrons, and the real one is emitted from the ini- 
tial state. The differential rate can be cast into the product 
of a leptonic and a hadronic tensor and the corresponding 
factorized phase space 



d(J = ^Lf,^Ht'''d<p2{pi,P2;Q,ki)d'Pn{Q;qi,-,q„)^ 



2tt 



(3) 



^ We disagree, however, with eq.(30) of the translated ver- 
sion [p3| which contains a missprint. 
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Fig. 1. Initial-state radiation in the annihilation process 
e+e~ 7+ hadrons at the Born level. 
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Fig. 2. One-loop corrections to initial-state radiation in e^e 
annihilation with the emission of a virtual photon. 



where d<l>n{Q] 9i, •, 9n) denotes the hadronic n-body phase 
space, including all the statistical factors coming from the 
hadronic final state. 

For an arbitrary hadronic final state, the matrix ele- 
ment for the diagrams in Fig. |l| can be written 



Mo 



where is the hadronic current and Af^ is the leptonic 
current in lowest order. Summing over the polarizations 
of the real photon, averaging over the polarizations of 
the initial e+e^ state, and using current conservation, 
Qi_iJ^ = 0, the leptonic tensor 



^0 



can be written in the following form: 

{Aiia/sf \ (2m^q^{l-q^f 2q^ 



q" 

8to2 



vl 



v\vl 



4<z2 



8m^\ P^iP^2+PiPf^ 



8m^ 



Aq^ 

J/12/2 



P^2P2 



(4) 
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with 



2ki ■ Pi 



m 



s 



-■ (5) 



It is symmetric under the exchange of the electron and the 
positron momenta. Expressing the bihnear products yi by 
the photon emission angle in the cms frame 



yi,2 



1 - 



(1 ^ 13 cost 



and rewriting the two-body phase space 

1 - g2 



327r2 



(6) 



it is evident that expression (^) contains several singular- 
ities: soft singularities for ^ 1 and coUinear singulari- 
ties for CO8 ^ ±1. The former are avoided by requiring 
a minimal photon energy. The latter are regulated by the 
electron mass. 

The physics of the hadronic system, whose description 
is model-dependent, enters only through the hadronic ten- 
sor 



(7) 



where the hadronic current has to be parametrized through 
form factors. For two charged pions in the final state, the 
current 



(8) 



where g7r+ E^^id g^- are the momenta of the 7r+ and tt^ 
respectively, is determined by only one function, the pion 
form factor [p^ . The hadronic current for four pions 
exhibits a more complicated structure and has been dis- 
cussed in [p|[25|,|26[|. 

At NLO, the leptonic tensor receives contributions both 
from one-loop corrections arising from the insertion of vir- 
tual photon lines in the tree diagrams of Fig. |^ and from 
the emission of an extra real photon from the initial state. 
In this paper, we consider only the emission of soft pho- 
tons. The implementation of these results in the program 
PHOKHARA and the discussion of their physical conse- 
quences will be considered in a separate work [ p7| . 

At NLO, the leptonic tensor has the following general 
forrnj^: 



2 r 



ai2 



an 



PlPl 



P1P2 



0.22 



P2P1 



(9) 



Terms proportional to are absent as a consequence of 
current conservation. The scalar coefficients aij and a-i 
allow the following expansion 



(0) 



a (1) 
— a) 



1-j 



a (1) 
a_i — — a_{ 

TT 



(10) 



The present convention differs from the one in Ref. by 
a factor 1/(1/11/2). 



The LO coefficients a|^'' can be directly read from eq. 



(0) 

^00 



"11 



"12 



2m2g2(l_g2^2 2q^ + yf+y^ 

yiy2 

^22 - On (yi 



8m^ 
If 



2 2 

yfyi 

4g2 



2/12/2 



y2j 



8m^ 
2/12/2 



(11) 



The imaginary antisymmetric piece proportional to a_i 
appears for the first time at NLO. The leptonic tensor 
therefore remains fully symmetric only at LO. 

As an alternative one can replace the Cartesian basis 
(eq. (H)) by a basis derived from the three circular polar- 
ization vectors of the virtual photon el and e±, defined 
through 



2 {q^'q•' - g^"' q^) 



■\/s g^(l — q^) 
2 {qf^q" - g^'^q^) 



kill J 



Vs 9^(1 - g2)2/3sin6l ' 

X ((2/1 - q^y2)piv - (2/2 - q^yi)p2v) 

^i'ripa'^lPlP2 I 



q^{l ~ q'^)(3 sin ( 



(12) 



where eii-qpa is the four-dimensional totally antisymmetric 
Levi-Civita tensor, £0123 = ^1- The leptonic tensor is thus 
given by 



i,j = L,±. (13) 



Only four of the scalar coefficients are independent 

Ol- aL+ , a-L = a+L = 0!+ , 
a = a++ , a |- = a-^-- . 

The trace of the leptonic tensor 

L'^'iq^q. - g^.q') = ^^^^^{aLL + 2a++) (14) 

is related to the cross section after angular averaging of 
the hadronic tensor. 

The relations between the components in the Cartesian 
and the circular basis read as follows: 



1 



All 



aL+ 



■[(2/2 - q^yifai 



4<72(l-q2)2 

+ (2/1 - '7^2/2)^022 + 2(yi - g^2/2)(2/2 - q^yi)a 



12 



4v/2q2(l - g2^ 
+ (1 + g^)(2/i - 2/2)012 - iira^i 
(3^ sin^ 9 



[(2/2 - g^2/i)aii - (2/1 - 9^2/2)022 



8 



-(an -I- a22 - 2ai2) 



a++ — — — aoo 



(15) 
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Conversely 



3 The NLO leptonic tensor 



an 



(l-g2)2 

(1 - g2)/3sin6' 

(1 -g2)2^2 sin^i 
022 = aii{yi ^ 2/2) , 
4v/V 



a_i = - 



Re(aL+) 



ai2 



an + a22 



4 04 



/32 sin^ 



7r(l - g2)/3sin( 



Im(aL+) 



(16) 



Combining the one-loop and the soft contribution we now 
arrive at the leptonic tensor in NLO. It will be convenient 
to split the coefficients a\^^ into a part that contributes 
at large angles and a part proportional to and m^, 
denoted by a^j''''' and a|j'™'' respectively: 



,(0) 



3 m 
2 (7^ 



log(4u;^)[l + log(m^)] 
2 ^2 " 



,(1,0) 



(l,m) 



(18) 



The scalar coefficients in the circular basis are given at 
LO by 



(0) _ 2g2/32 gjj^. 



yij/2 



^(?/i - y2)/3sin( 

V2{1 - g2)yiy2 



2m2(l - g2)i 



(0) ^^sin^/ 2 , m2(l-g2)^ 
all = I -q H 



yi2/2 V yiy^ 

(0) _ (0) , 2q^ + yl + yj 2m^q^{l-q^f 



2/12/2 



9 9 

2/12/2 



(17) 



The factor proportional to the LO coefficients a-^"* con- 
tains the usual soft and coUinear logarithms. The expres- 
sions are particularly compact in the circular basis. For 

completeness we also repeat the results for a-j'"'', which 
can be found in albeit in the Cartesian basis: 



,(1,0) 



(l-g2)2 



yi 



1 - yi 

„2 



21og(g2) + 2i(yi) 



1 - y2 (1 - 2/2)^ 



log(^) + [2/1 



2/2 J 



(19) 



The one-loop matrix elements (Fig. ||) contribute to 
the leptonic tensor through their interference with the low- 
est order diagrams (Fig. Q). They contain ultraviolet (UV) 
and infrared (IR) divergences, which are regularized using 
dimensional regularization m D = A — 2e dimensions. The 
UV divergences are renormalized in the on-shell scheme. 
The IR divergences are cancelled by adding the contribu- 
tion of an extra soft photon emitted from the initial state 
and integrated in the phase space up to an energy cut- 
off < Wy^ far below ^/s. The result, which is finite, 
depends, however, on this soft photon cutoff. Only the 
contribution from hard photons with energy Ey > 
would cancel this dependence. 

The algebraic manipulations have been carried out with 
the help of the FeynCalc Mathematica package |2^. Us- 
ing standard techniques [29[| , it automatically reduces the 
evaluation of the one-loop contribution to the calculation 
of a few scalar one-loop integrals and performs the Dirac 
algebra. 

Since we consider the small angular region, mass terms 
proportional to y~'^ and even y^^ arise. Terms propor- 
tional to m? and even must be kept, if they are mul- 
tipHed by y~'^ and y~^ respectively. In the expansion of 
the one-loop integrals, functions that depend on the ratio 
m?' /yi cannot be expanded, in contrast to functions of 
or yi separately. To arrive at a systematic approach we 



therefore make the replacements m 



Am 



2/* 



A2/j, 



perform the expansion for small A up to the appropriate 
order, and set A = 1 at the end. 



,(1,0) _ 



,(1,0) 



/3sin( 



2/1(1 + g')' 
2) i(l-yi)(l-y2) 



log(g') 



2g^(2 + g^) 

2/2 

2(1 + q^{l + q^Y 



2(l-g^) --^ 
2/1 2/2 



(1-gT 

2/2 
2 



2q^ 
1 - 



1-2/2 (1 - 2/2) 
l + 2g2 



log(^) 
q^ 



2/1 



yl 



1 



2/1 



(l-2/i)^ 



i(2/i) 
- [2/1 



mq 



2 r 



2/12/2 
2/2] 



21o g(l-2/i) 
2/1 



(20) 



1/1-272 



1 



2/1 



1 + g^ 
{\-q^Y\\-y^ 



-2 + 21og(g^) + 



1 



l-g2 
2 

l-g2 



1 + g^ 
\-q^\(\-y^f 
(l + g')2/2 , 2 
(1 - 9^)2/1 2/2 



- 1 



g 



1-2/2 

[2/1 ^ 2/2] 



log — 

g 



(21) 



^ The result for the imaginary part of aL+ differs from a_i 
in the original version of ]lq]. 
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[l-q^f q\l-q^) 



Vi 



log(g2 



V2 



X log — , 

q 



(1 



1 - y2 2(1 - 



Hvi) + [yi ^ 2/2] 



(22) 



that the mass-suppressed terms are all real. The functions 
n{yi, z) and N{yi) are defined through 



n{yt,z) 



yii-m? - yi) 



1 + z log(4) 



[m? - yi)'- 



iog(4) 



(27) 



with 



L{y{) ^ Li2{^) - U2{1 



l)+log(g2 + y01og(^) 
q'^ 



where Li2 is the Spence or dilogarithm function. The re- 
sult in the Cartesian basis has been given in Ref. . The 
terms proportional to powers of are given by 



(l.m) 



= 



(23) 



,(1,™) 



1(^13 sin J 4m^ 



vl 



log(g2)log(^) 



+ ALi2{l~q') + -[Li2{l 

2 V m 



n{yi 



2(1 -3g2 



^1 _ _ 

2> g 

2m^N{yi) 



yi(m2(l - g2) „ j^^) 



[yi <-> y2] L 



(24) 



(l.m) _ _toV 
^++ 2 



Iog(q^)log(^)+4iz2(l-g2) 
m q 



Lt2(l ^) - — 

m- D 



i-iog(^) 



-("^^2(1-^)-- 

[yi ^ ^2] , 



(25) 



7V(2/,) = log(g^) log(^) + Li2{l-q') 
+ Li2{l 2) - IT ■ 

TO^ D 



(28) 



The apparent singularity of n inside the limits of phase 
space is compensated by the zero in the numerator. For 
the numerical evaluation in the region yi close to we 
use 



n{yi,z)\ 



1 + z log(^) 



n— 



(29) 



For the conversion from the circular to the Cartesian 
basis, and to ensure finite results in the limit smd — )■ 0, it 

is important that a l+ and a.] vanish at ~ sin 9 and sin^ 9 

respectively. This corresponds to the requirement that the 
factors in curly brackets do not diverge for small sin0, i.e. 
in the limit {m?{l ~ q^) — yi) ^ 0. This is guaranteed by 
the behaviour of N{yi) for yi m^(l — q^): 



m'^N{yi) 



m2(l - g2) _ y. 



log(l - q^) \og{q^ 



yi—^m'^ (1— q2 



1 



(30) 



and 



(l,m) 



sin 9 r Am 



yl 



log(g2) log( 



yt . 

m'^q^' 



+ 



2(u2il-^)-^- 
\ 6 

2TO^(l-g2) 
?/i(m2(l - g2) _ y^) 



■log( 



q 
yi 



The results for a^j''"'' in the Cartesian basis are listed in 
iLi2{l-q^) Appendix 0. 

We note that the imaginary part of i^^, which is 
present in the coefficients a^^ or a_i only, is of interest 
for those cases where the hadronic current receives contri- 
butions from different amplitudes with non-trivial relative 
phases. This is possible, e.g. for final states with three or 
more mesons or for pp production. 



— —n(yi, 



1-^2 

log(g^) 
l-g2 



-1- 



3 - ^2 

— I 

1 — g2 m2(l — q'^) — yi 

[yi ^ 2/2 



N{yi) 



(26) 4 Tests of the result 



The coefficient a^^ is antisymmetric with respect to the 
exchange [yi y2]^ while all the others are symmet- 
ric. Only a^^_j''^ contributes to the imaginary part. Notice 



After summation over the polarizations of the virtual pho- 
ton the differential rate is given by the properly contracted 
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leptonic tensor: 

,v2 



2 L^^^iq^qy - g^^uq^) = aLL + 2a++ 



= -3 aoo + (^(1 - yi f - "^^) 0,11 
+ (^(l-y2)'-m2)a22 
+ (7^(l-yi)(l-y2)-(l-2m2))ai 



^2q 



log(4w^)[l + log(m^)] 



3, „ TT^ 



2(1 



.log(g2 



1 - 



2/12/2 



1 



X log(^) - 2^ ' Afi 

q^ 2(1-9^) 4(1 



[yi <-> 2/2] ) - 2a^^+"^ 



2/1 - 2^2 

2(1 -2;2) 4(l-y2)2 
1 1-2/2 



42/1 



(31) 



where 



(47ra/s)2 
2 



Lo^iqf^q^ - q^9i,u) 



2/12/2 



2g' + 2/? + 2/2 



2mV(i-g')' 



2/12/2 



(32) 



Lg" being the leptonic tensor at LO. In Lq only the rel- 
evant terms in the limit m? have been kept. Large 
angle terms and mass corrections are in agreement with 
Berends et al. [^|2^. After proper analytic continuation 
and using the substitutions 



t 



-2/1 + iV 



-2/2 + iv 



s + iri , (33) 



the leptonic tensor in eq. (^ is closely related^ to the ten- 
sor T^i/, which describes "virtual Compton scattering" and 
has been calculated by Kuraev et al. p2| . After performing 
this analytical continuation, the results are in agreement 
(However, the result for T12 printed in contains a typo 
and would give rise to a discrepancy.). 



5 Conclusions 



Compact analytical formulae have been obtained for the 



virtual corrections to the process e"^ 



77*, which are 



valid for photon emission under both large and small an- 
gles. After proper analytic continuation the results are in 
agreement with those obtained in for the reaction 
"/* + e~ 7 + e~. In polarization averaged form they 
are in agreement with those for e+e^ — > 7Z from ||20|,pT|. 

We thank N.P. Merenkov for drawing our attention to this 
reference. 



Acknowledgements 



We would like to thank G. Cataldi, A. Denig, S. Di Falco, 
W. Kluge, S. Miiller, G. Venanzoni and B. Valeriani for 
reminding us constantly of the importance of this work for 
the experimental analysis and H. Czyz and N.P. Merenkov 
for very useful discussions. Work supported in part by 
BMBF under grant number 05HT9VKB0 and E.U. EURO- 
DA^>NE network TMR project FMRX-CT98-0169. 



A The leptonic tensor in the Cartesian basis 



For the convenience of the reader we list the mass-suppressed 
terms ajj'™' also in the Cartesian basis. The large-angle 
contributions have been given in . The component pro- 
portional to (7^1/ reads 



(l,m) 



(34) 



see eq. (25). The coefficient does not contribute at 

this order to aoo- The remaining components are given by 



(l,m) 
Hi 



q^ ( 4m^ 



l-log(^) 



Ui \ 



"(2/1, 1) 



2m2g2 



yi(m2(l - g2) _ 
2 

III. 

1 



m 



to2(1 _ q2-j _ 
1 r4TO2(l-g2) 



1 1 ^ 2/1 ^ . log(g^) 
— log — + -j 2 



N{yi) 



log(g2)log(^) 



4L*2(1-9') + 2(Lz2(1-^)-^ 

i-iog(4) + (i + -)^^2(i-4) 

77?/ \ 2/2 / \ ^ 

2— '^(^2, —4 ) 



4771^9^ 


2/1 













2m? 



y2(m2(l - g2) _ ^3) 

9 

m 

3 



7772(1 - q^) - y2 



(l,m) (li^Ti) 



1 - 2g4 
4log(4) + ^ 

qz 1 — 

N{y2) 



*22 



Oil '{yi ^ 2/2) 



(35) 



(36) 
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,(1,™) 

*12 



Arrr 



l-g2 

' 1 2 

1 m 

2 ^ y7 



1 - log(^) 



are obtained: 



yi 



2 \ 1 



COl = 



2m2 



C02 = 



i«2(l J-) - -r 



zi-21og(/3)- 



2^2 



log(c) 



log(c) 



2Lz2(c)- — +i7r(Z\-21og(/3) 



777,2(1 — q^) — yi 



N{yi) 



[yi ^ 2/2] 



C04 = 



(37) 



1-^2 

-2 



i^„i^+,,log(^) 
2 2 ^^Cg 



^3y^ 



Z\ + 21og(-) log(c5)+log^(c) 



and 



+ 2i72(l - -) + 2X^2(1 -CCg)- Ll2{l - C^) 



aLT^ - 



(38) 



77r Z\ + 21og( ' 

cVi' 



B Scalar one-loop integrals 



A few two-, three-, and four-point scalar one-loop inte- 
grals enter our calculation. Expression for the two-point 
scalar integrals are simple and well known. The notation 
from [^8|, where the corresponding results valid for large 
photon angles (7772 <^ 1, q^,yi) can be found, is used in the 
following. The general three-point scalar one-loop integral 
is defined by 

U){Pa, [Pa-Pbj ,Ph,m^, 171^,171^) = -7l67r ^ 

d^k 1 
{2tt)D [fc2 _ ^2] [(fc - p,)2 _ ^2] [(fc _ p,)2 _ ^2] ' 

(39) 



Four different three-point scalar one-loop integrals are needed 

COl = Co((p,-A:i)2, 0,777^,0, 7772,777^) , 
C02 = Co(777g, s, 777g, 0, ml, ml) , 
COS = Co{ip^ - ki f,Q^,ml,0,ml,ml) 
COA ^ CoiQ^ , s,0,ml,ml,ml) , 



with 



and 



A = 



(42) 



(43) 



P = y/l - 47772 



'1 



4r7i2 



1-/3 
1 + 



(44) 



Our expression for the C03 function is rather cumbersome: 
r /(i-zi)z2~ 



C03 = 



g2(zi - Z2) 



log(^)log 
9 



(1 - Z2)zi 



+ < Li2{ — ) - log((z3 - Zi)(z4 - Zi)) l0g(l - — ) 



+ Li2{ 



Zl 

Zl' 



)-Lt2{- 



'Zl 



Li2{- —) - Li2{ 



where 



i = 1,2, and one scalar box 



DO 



2,,4-D 



(2^) 



D 



(40) 



(41) 



21.2 



Z3A 



Z4 - Zl 



q^ + yx 



Z3 - Zl 

Zl - 
■ Z4 - Zl ' 



[Zl ^ Z2\ 



(45) 



2g2 



i± Wi- 



4(7772 _ ij^^q2 



(<?' + y^V 



l±^ I 



4(7772 _ jj^-j 



being 



X 



k^[{k + p^)^ - ml][{k + p^ ~ kiy - ml][{k ~ pj)^ ~ m^] ' /to(C03) 



(^^3 - Zi)(z4 - Z2) 
g^(zi - Z2) V(24 - Zi)(z3 - Z2) 



log 



(46) 



(47) 



with j y^i. 

its imaginary part. As for the other scalar one- loop func- 
The following simple expressions are used, from where tions, this expression is expanded up to terms of order 7772 
the limits (7772 ^ l,q'^,yi) or (7772 <c 1,?^ but 7772 ~ yi) and 771'*. 
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